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Abstract —  In  this  work,  we  perform  a  comparative  study  of  flow  models  and  a  sensitivity  analysis  con¬ 
cerning  some  parameters,  considering  the  effects  of  slipping  and  gas  adsorption  in  shale  gas  reservoirs, 
using  numerical  simulation.  We  use  the  Finite  Difference  Method,  a  linearization  based  on  the  Picard 
method,  and  the  Conjugate  Gradient  method  to  obtain  the  reservoir  and  production  well  pressures.  The 
results  obtained  demonstrate  the  importance  of  the  effects  of  slipping  and  adsorption  in  the  variation  of 
pressure  in  a  single-phase  isothermal  gas  flow.  Mainly  because  some  parameters  in  the  flow  models, 
depending  on  their  values,  can  act  favoring  (or  disfavoring)  the  pressure  drop. 

Keywords  —  Adsorption.  Finite  Difference  Method.  Naturai  gas.  Reservoir  simuiation.  Siipping 
fiow. 


I.  INTRODUCTION 


face  conditions),  for  exampie,  when  composed  essen- 
tiaiiy  of  methane.  Generaiiy,  naturai  gas  (depend¬ 
ing  on  the  reservoir  and  its  characteristics)  contains 
methane,  ethane,  propane,  butane,  pentane,  and  hex¬ 
ane,  but  we  can  aiso  find  other  heavier  components  in 
it.  The  non-hydrocarbon  parceis,  considered  as  impu¬ 
rities,  inciude  carbon  dioxide,  hydrogen  suifide,  and 
nitrogen.  The  primary  component  is  aiways  methane, 
which  may  represent  70  to  98%  of  the  totai,  foiiowed 
by  ethane  (which  may  reach  10%). 


For  decades,  the  oii  industry  has  foiiowed  a  stan¬ 
dard  strategy  in  the  context  of  hydrocarbon  expio- 
ration/expioitation:  the  search  for  rock  formations  of 
interest,  the  iocation  of  the  reservoir  rocks,  the  identi¬ 
fication  of  the  mechanisms  of  imprisonment,  and  weii 
driiiing.  These  procedures  had  been  used  by  geoi- 
ogists,  geophysicists,  and  engineers,  aiming  mainiy 
at  the  production  of  hydrocarbons  via  verticai  weiis  in 
reservoirs  formed,  for  exampie,  by  sandstones  or  car¬ 
bonates. 


Naturai  gas  is  the  cieanest  and  richest  source 
of  hydrogen  among  hydrocarbons  and  has  high  effi¬ 
ciency  in  energy  conversion.  The  current  trend  is  for 
naturai  gas  to  become  one  of  the  most  important  fu- 
eis  in  the  giobai  economy  [12]  due  to  economic  and 
environmentai  reasons  [29,  31].  See  Fig.  1  for  the 
breakdown  of  the  internai  energy  suppiy  in  Brazii  for 
the  year  2019  [13]. 


However,  in  the  past  two  decades,  new  concepts 
of  expioitation  have  emerged.  Rock  formations  that 
were  previousiy  iess  important,  to  the  point  of  not 
being  considered  economicaiiy  viabie,  are  currentiy 
seen  as  potentiai  reservoirs.  Among  unconventionai 
sources,  shaie  gas  reservoirs  are  gaining  more  and 
more  reievance  on  the  worid  stage  [8]. 


Naturai  gas  is  a  homogeneous  fluid  with  iow  vis¬ 
cosity  and  density.  In  a  reservoir,  naturai  gas  can 
be  ciassified  as  dry  (it  does  not  contain  moiecuies 
heavy  enough  to  form  a  iiquid  hydrocarbon  under  sur¬ 


1.1.  Naturai  gas 


In  terms  of  its  use,  naturai  gas  can  provide  heat¬ 
ing,  generate  eiectricity,  and  we  can  use  it  as  fuei  for 
motor  vehicies  or  in  the  production  of  hydrogen  [12]. 
According  to  Economides  and  Wood  [12],  the  worid 
reserves  of  naturai  gas  have  been  increasing  at  an  an- 
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nual  rate  of  about  5%  since  the  seventies.  The  num¬ 
ber  of  countries  with  known  naturai  gas  reserves  has 
grown  from  40,  in  1960,  to  85  in  2005. 

1.2.  Shaie  gas 

Non-conventionai  reservoirs  differ  from  conven- 
tionai  ones,  for  exampie,  in  that  they  have  a  very 
iow  permeabiiity.  Among  the  non-conventionai  reser¬ 
voirs,  we  can  mention  the  tight-gas  sands  (iow  per¬ 
meabiiity  reservoirs),  the  coai  bed  methane  (in  which 
there  is  gas  adhered  to  coai  veins),  and  the  shaie  gas 
(gas  in  shaie-iike  rocks).  The  iatter,  the  focus  of  this 
work,  typicaiiy  have  extremeiy  iow  permeabiiities,  in 
the  range  of  10~®  to  10~®  Darcy,  iow  porosity,  gas 
adsorbed  on  the  soiid,  and  the  gas  siips  on  the  pore 
surfaces  of  the  reservoir.  Its  formations  are  composed 
of  sedimentary  rocks,  consisting  of  ciay  minerais  such 
as  iiiite,  smectite,  and  kaoiinite  [8]. 

Energy  Matrix  2020 
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the  pores,  dissoived  in  iiquids,  or  adsorbed  on  soiids. 
We  find  the  gas  adsorbed  on  carbon-rich  compounds, 
such  as  kerogen.  According  to  Aii  [1],  Wang  [27]  and 
Berawaia  [3],  gas  in  an  adsorbed  state  represents  a 
fraction  ranging  from  20  to  80  %  of  the  totai  gas  re¬ 
serves  in  shaies.  As  the  pressure  decreases,  as  a 
resuit  of  continuous  production,  the  gas  adsorbed  de¬ 
taches  from  the  soiid  and  becomes  part  of  the  free 
phase,  contributing  to  flow  and  production.  Therefore, 
by  negiecting  the  effects  of  adsorption,  we  can  under¬ 
estimate  the  potentiai  of  a  shaie  reservoir  [30]. 

About  two  decades  ago,  we  considered  gas  re¬ 
covery  in  shaie  reservoirs  as  economicaiiy  and  tech- 
nicaiiy  nonviabie.  However,  with  the  emergence  of 
new  technoiogies,  especiaiiy  the  driiiing  of  horizontai 
weiis  and  hydrauiic  fracturing,  their  expioitation  has 
become  possibie  and  economicaiiy  viabie.  In  2015, 
about  three-quarters  of  natural  gas  production  in  the 
United  States  originated  from  these  reservoirs,  with 
the  expectation  of  continued  growth  over  the  next  two 
decades  [28].  The  great  success  in  North  America 
has  led  to  unconventional  sources  gaining  more  at¬ 
tention  on  the  world  stage,  including  the  countries  of 
South  America,  such  as  Argentina  and  Brazil. 

II.  NATURAL  GAS  FLOW  IN  POROUS 
MEDIA 

Typically,  the  classical  Darcy  law  [7]  is  used  in  en¬ 
gineering  to  describe  low-speed  flows  in  porous  me¬ 
dia.  For  non-Darcy  flows,  it  is  modified  and  expressed 
in  the  form: 

V  = -— {Vp- pgVD) ,  (1) 

where  v  is  the  surface  velocity  of  the  fluid,  ka  is  the 
apparent  permeability  tensor,  p  is  the  viscosity,  p  is 
the  pressure,  p  is  the  density,  g  is  the  magnitude  of 
the  acceleration  due  to  gravity,  and  D  is  the  depth. 
Among  the  non-Darcy  effects,  which  we  can  incorpo¬ 
rate  into  the  apparent  permeability,  we  have: 

•  the  inertial  and  turbulent  effects, 

•  the  non-Newtonian  flow. 


Fig.  1 :  Supply  by  energy  source  in  Brazil  in  2019. 


•  the  slipping  and  adsorption  of  the  gas  on  the 
pore  surface. 


In  them,  the  gas  is  present  as  a  free  phase  in  in  addition  to  others  that  also  cause  non-linearity  [2]. 
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2.1 .  The  gas  slippage  phenomenon 

Florence  et  al.  [15]  reported  that  the  phenomenon 
of  gas  slipping  occurs  when  the  mean  free  path  of  the 
fluid  molecules  is  of  the  order  of  the  characteristic  hy¬ 
draulic  radius  of  the  pores.  We  know  that  gas  flows 
in  the  porous  medium  differently  from  a  liquid  for  two 
reasons:  its  high  compressibility  and  the  Klinkenberg 
effect  [2].  We  use  Klinkenberg’s  correction  to  take  into 
account  the  effects  of  gas  slippage,  which  can  be  veri¬ 
fied  when  carrying  out  permeability  measurements  on 
rock  samples  [14] 


where  b  is  the  Klinkenberg  parameter  and  k  the  abso¬ 
lute  permeability  tensor  [19].  For  example,  gas  slip¬ 
page  can  occur  in  shale  gas  and  tight  gas  sand  reser¬ 
voirs. 

As  an  alternative  to  Klinkenberg’s  correction,  we 
can  also  account  for  this  effect  through  a  more  general 
expression  for  the  determination  of  apparent  perme¬ 
ability,  which  we  can  calculate  as  a  function  of  Knud- 
sen  number  [18], 


of  the  porous  medium,  Z  is  the  compressibility  fac¬ 
tor,  T  is  the  temperature,  R  is  the  universal  constant 
for  gases,  M  is  the  molecular  mass  of  the  gas,  and 
7  =  M/Mair  is  the  relative  density  of  the  gas  and 
Mair  the  molecular  mass  of  air  [8].  The  slip  regime 
exists  for  the  10~^  <  Kn  <  0.1  range  [4]. 

2.2.  The  adsorption  phenomenon 

Adsorption  occurs  when  the  gas  molecules  attach 
to  the  solid  surface  of  the  pores,  and  the  volume  of 
gas  adsorbed  depends  on  the  pressure  of  the  gas 
phase.  Here,  we  consider  the  Langmuir  isotherm 
model: 


Tarfs  — 


pVl 


U) 


P  +  Pl 

where  Vads  is  the  volume  of  adsorbed  gas,  14  is  the 
maximum  volume  adsorbed  (Langmuir  volume),  and 
Pl  is  the  pressure  corresponding  to  half  of  the  max¬ 
imum  adsorption  capacity  (Langmuir  pressure).  Ac¬ 
cording  to  Jiang  and  Younis  [18],  we  can  also  incor¬ 
porate  the  effects  of  adsorption  into  the  calculation  of 
apparent  permeability 


/  AKv  \ 

ka  = /(A:n)k  =  (1 -h  offcATn)  f  1 -h  j  k,  (3) 

where  k  is  the  absolute  permeability  tensor,  and 
is  the  rarefaction  parameter  [18],  whereas  for  the  slip 
regime,  au  =  0.  The  Knudsen  number,  Kn,  is  defined 
by 


(  L.  \ 

^  dm 

PL 

Rh 

\  PL  / 

where  is  the  diameter  of  the  gas  molecule  adhered 
to  the  surface.  Therefore,  we  can  combine  the  contri¬ 
butions  due  to  gas  slippage  and  adsorption  in  deter¬ 
mining  the  apparent  permeability  in  the  form: 


Kn 


A 


(4) 


k,  =  f{Kn')\^a 


(9) 


where  A  is  the  average  free  path  of  the  molecules  and  where 
Rh  the  characteristic  hydraulic  radius  [15] 

/(AV)  =  (l  +  atAV)(^i±^)  (10) 

and  the  modified  Knudsen  number  is  set  to 


Rh  =  2V^J  - 


and 


p  I  ttZRT 
^  p\  2M 


(6) 


where  t  is  the  tortuosity  of  the  porous  medium,  k  is 
the  geometric  mean  representing  the  absolute  perme¬ 
ability  of  the  porous  medium,  4>  the  effective  porosity 


Kn'  = 


A 

^eff 


where 


/  ^ 
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is  the  effective  hydrauiic  radius. 

2.3.  Governing  equation 

According  to  Li  et  ai.  [21],  the  mass  conservation 
equation,  incorporating  the  effects  of  adsorption,  is 
given  by 

^  ( Psc^'\  d  f  PsPscVads\  yy  /  PscV\ 

B  \~ir) 


where  Qsc  is  a  source  term,  ps  is  the  density  of  the 
rock,  B  =  psc/p  is  the  voiume  formation  factor,  14 
is  the  totai  voiume  of  the  voiume  controi  (rock  pius 
pores),  and  the  subscript  sc  indicates  the  standard 
conditions, 

0  =  ^0  [l  +  C0(p-/)]  ,  (14) 


being  that  we  used  the  term  source  to  represent  the 
weii  flow  through  an  internai  boundary  condition  for 
the  weii-reservoir  coupiing: 

Qsc  -  Jw  (/P  Pwf')  ?  (1  8) 

where  J^,  is  the  productivity  index,  and  p^f  is  the  pres¬ 
sure  in  the  weii  [14].  In  this  work,  we  empioy  a  pre¬ 
scribed  production  flow  condition. 

We  know  that  Eq.  (15)  is  a  partiai  noniinear 
paraboiic  differentiai  equation  that  we  use  to  deter¬ 
mine  the  gas  pressure.  As  an  initiai  condition,  we 
consider 

p{x,y,t  =  0)  =  pini{x,y),  (19) 

where  p™  is  the  initiai  pressure  before  the  reservoir 
is  disturbed  by  production/injection. 

On  the  other  hand,  the  externai  boundary  condi¬ 
tions  are  those  of  nuii  flow  at  the  borders 


where  the  superscript  0  indicates  the  reference  con¬ 
ditions,  and  C0  is  the  compressibiiity  coefficient  of  the 
rock.  We  consider  here  smaii  and  constant  rock  com¬ 
pressibiiity. 

Substituting  Eq.  (1)  in  Eq.  (13)  and  considering 
that:  psc  and  ps  are  constant,  we  negiect  the  gravi- 
tationai  effects,  the  flow  is  two-dimensionai,  and  the 
apparent  permeabiiity  tensor  is  diagonai. 


/  kgy  dp\ 

dx  \  pB  dx  )  dy  \  pB  dy  ) 


Vb 


(p-Pwf) 


V  =  0,Ly 


=  0, 


(20) 


where  L^.  and  Ly  are  the  respective  iengths  of  the 
reservoir  in  the  x-  and  y-  directions. 


III.  NUMERICAL  METHODOLOGY 

In  this  work,  we  use  the  Finite  Difference  Method 
(FDM)  and  a  computationai  mesh  of  centered  biocks 
[2,  14],  together  with  a  weii-reservoir  coupiing  tech¬ 
nique  [1 0],  to  determine  the  pressures  in  the  reservoir 
and  the  producing  weii. 


=  (r;+n)|  (15) 

where  for  a  non-isothermai  flow  and  considering  the 
previous  expressions  of  14ds  and  B 


r' 


and 


(16) 


r;  =  Ps 


1  d^ads  I  ri 

B  dp  ^ 


(17) 


3.1 .  Discretization  of  the  governing  equation 

We  provide  a  schematic  representation  of  a  dis¬ 
cretized  two-dimensionai  domain  is  provided  consid¬ 
ering  a  singie  ceii  in  the  direction  of  the  z-axis  (high- 
iighted  in  Fig.  2).  We  obtain  the  numericai  soiution  at 
the  nodes  of  the  computationai  mesh,  iocated  in  the 
centers  of  the  ceiis.  Ux  and  Uy  represent  the  number 
of  ceiis  in  the  x-  and  y-  directions,  respectiveiy.  The 
integer  indexes  i  and  j  represent  the  numbering  of  the 
ceii  nodes  in  the  respective  x-  and  y-  directions,  and 
the  fractionai  indexes  i  ±  1/2  and  j  ±  1/2  their  inter¬ 
faces. 

Therefore,  by  writing  the  governing  equation  as 
being  evaiuated  in  ceii  node  i,j,  and  at  time  ievei  n+1: 
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1  n+1 


where,  as  +  =  dxdyL^,  we  used 


(21) 


n+1 


(25) 


Spatial  derivatives,  evaluated  on  the  faces  of  cells, 

(Fb)ij  =  {^x^y)ij  Lz  (22)  are  also  discretized  by  centered  differences, 

and  we  introduced  the  variables 


T'  = 


A  h 


and 


(23) 


A  h 

™/  _ 
y  ~  yB 


(24) 


with  +  =  AyLz:  and  +  =  AxL^.  Axij  and  Ay^j 

are,  respectively,  the  mesh  spacing  in  the  x-  and  y- 
directions  in  the  cell  and  the  length  of  the 
rock  formation  in  the  z-direction. 


>j 


Fig.  2:  Example  of  some  cells  of  the  discretized 
domain. 

Then,  we  employ  a  centered  difference  type 
scheme  to  approximate  the  spatial  derivative  in  the 
^-direction  [2,  14], 


I  ^  \  ^  P»+l,j  -Pjj 


and 


(26) 


d"+1 


n+1 


Ax, 


(27) 


where  Axi±i/2j  is  the  distance  between  the  centers 
of  cells  (i,j)  and  (i  ±  1,  j).  We  can  obtain  the  approx¬ 
imations  for  derivatives  in  the  y-direction  in  an  analo¬ 
gous  way. 

We  also  introduced  transmissibilities  in  the  x-  and 
y-  directions: 


(28) 


and 


■J'n-l-l  _  f  d^yka,v\ 


(29) 


where  a  harmonic  mean  is  used  to  determine  the  val¬ 
ues  of  areas  and  permeabilities  in  positions  {i±l/2,j) 
and  {i,j  ±  1/2),  from  the  values  known  in  neighboring 
nodes.  With  regard  to  the  properties  of  the  fluid,  we 
applied  an  arithmetic  mean  [14]. 

Thus,  using  an  Euler  approximation  for  the  time 
derivative  and  a  fully  implicit  formulation  in  time,  it  is 
possible  to  obtain  the  final  discretized  form  of  Eq.  (21 ), 
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+  (qsc) 
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(30) 


where  n  indicates  the  time  ievei  at  which  we  know  the 
pressure,  and  we  used  a  conservative  expansion  for 
the  accumuiation  term  [14]: 


3.3.  Linearization 

Equation  (30)  forms  a  system  of  noniinear  aige- 
braic  equations  [17].  From  Eq.  (33),  for  the  particuiar 
case  of  the  two-dimensionai  probiem,  we  can  isoiate 
the  pressure  of  the  verticai  weii  once  we  know  the  flow 
and  the  pressure  of  the  reservoir  at  the  end  of  each 
time  step. 

To  achieve  the  iinearization  of  Eq.  (30),  we 
adopted  the  weii-known  method  of  Picard  [22].  There¬ 
fore,  after  its  iinearization,  we  get  the  foiiowing  equa¬ 
tion 
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3.2.  Weii-reservoir  coupiing 

We  proceed  with  the  process  of  discretization  and 
the  source  term  (gsc)"^^  takes  the  discrete  form 


/  \n+i  /  j  \n-\-i  n+1  /  \Ti 

(<?sc),  P,J  -  {Pwf)i 


\n+l 


(33) 


and  the  productivity  index  is  expressed  by  [24] 


(+) 
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'i'J 
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27rZ/2  yJka^xPa.y 


Bfi  In 


(34) 


where  is  the  radius  of  the  weii  and  the  equivaient 
radius,  r^q,  is  given  by  [24] 


and  we  aiready  rewrote  the  equivaient  radius  for  its 
use  with  the  apparent  permeabiiity. 
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(36) 


and  transmissibiiities,  Tp,  Tg,  and  Qsc  are  evaiuated  at 
the  iterative  ievei  v. 

Then,  we  utiiize  the  vaiues  determined  at  the  iter¬ 
ative  ievei  V  to  caicuiate  the  pressure  at  i;  -i-  l,n  -i- 1. 
Therefore,  we  obtain  the  pressure  from  two  inter- 
nai  and  externai  iterative  processes  for  a  given  time 
step  [8]. 

We  choose  the  method  of  the  Conjugate  Gradient 
(CG)  to  soive  the  aigebraic  system  of  Eq.  (36)  [1 6,  25]. 

We  can  find  in  Fig.  3  the  flowchart  of  the  numerical 
scheme  of  the  resolution  algorithm  for  a  single  time 
step  calculation.  To  verify  whether  we  achieve  the 
convergence  of  the  internal  and  external  iterations,  we 
use  tolerance  values  equal  to  toh  and  toh,  respec¬ 
tively. 


IV.  RESULTS 

In  this  section,  we  present  our  numerical  results. 
Initially,  we  introduce  the  standard  case  and  some  in¬ 
formation  about  the  data  used  in  the  simulations,  in- 
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_ JL _ 

End  of  time  step 


Fig.  3:  Flowchart  for  a  single  time  step. 
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eluding  rock  and  fluid  properties,  simulation  parame¬ 
ters,  and  geometric  data. 

Next,  we  do  a  numerical  verification  study,  includ¬ 
ing  the  refinement  of  the  computational  mesh.  Finally, 
we  carry  out  a  sensitivity  study  of  parameters  and 
compare  the  different  models  employed.  Here,  the 
standard  model  is  that  incorporating  the  effects  of  gas 
slippage  and  adsorption  in  the  modified  Darcy’s  law, 
as  discussed  in  Jiang  and  Younis  [18]. 

We  show  the  results  using  specialized  plots,  com¬ 
monly  employed  in  the  area  of  well  pressure  test  anal¬ 
ysis  [5]. 

We  consider  six  models  according  to  how  we  cal¬ 
culate  the  apparent  permeability  values  (ka)  in  Eq.  (1 ) 
and  the  coefficient  in  Eq.  (15): 

•  Model  1  [14],  which  represents  the  classic  case 
of  Darcy’s  law,  without  the  effects  of  slippage 
and  adsorption; 

•  Model  2  [15],  the  apparent  permeability  is  calcu¬ 
lated  as  a  function  of  the  Knudsen  number  (slip¬ 
page  effect); 

•  Model  3  [3],  without  slippage  effect  and  with  an 
accumulation  term  for  the  inclusion  of  adsorp¬ 
tion; 

•  Model  4  [21],  apparent  permeability  corrected 
as  a  function  of  the  Knudsen  number  (slippage 
effect)  and  insertion  of  adsorption  through  the 
source  term  in  Eq.  (1 5); 

•  Model  5,  we  modify  the  apparent  permeability 
taking  into  account  adsorption  and,  also,  we  in¬ 
troduce  its  effects  via  the  accumulation  term  in 
Eq.  (15); 

•  Model  6  [1 8],  we  change  the  value  of  to  take 
into  account  the  combined  effects  of  slippage 
and  adsorption,  and  we  also  include  the  term 
Fs  7^  0  in  Eq.  (15). 

We  modify  the  initial  time  step  (Ai™)  according  to 
the  growth  rate  (SAt)  to  obtain  the  next  time  step.  This 
procedure  is  interrupted  when  we  reach  the  specified 
final  time  step  (Atmax)-  We  employ  this  strategy  to 
achieve  greater  accuracy  in  the  initial  stages  of  pro¬ 
duction. 

We  can  find  the  default  values  used  in  the  simula¬ 
tions  in  Table  1 .  They  are  defined  based  on  the  works 


of  [21],  [18]  and  [9].  We  adopt  the  same  tolerance 
value  for  the  internal  (tok)  and  external  (tol2)  iterative 
procedures. 

Table  1 :  Parameters  for  the  default  case. 


Parameter 

Value 

Unit 

1.0x10-® 

psi-i 

dm 

2.3x10-1° 

ft 

5.0x10-° 

Darcy 

II 

4,000 

ft 

L, 

40.0 

ft 

Lwf 

40.0 

ft 

Tlx  —  '^y 

321 

- 

Psc 

14.65 

psi 

PL 

1,100 

psi 

Pini  =  P° 

6,000 

psi 

Q  sc 

-5.0x10^ 

scf/day 

^  w 

0.25 

ft 

R 

10.73 

ft°psi/R  Ibm-mol 

imax 

365 

day 

toll  =  tol2 

1.0x10-4 

psi 

T 

609.67 

R 

Tsc 

519.67 

R 

Vl 

0.0005 

ft°/lbm 

1 

0.6 

psi 

^At 

1.05 

- 

^^ini 

0.01 

day 

^^max 

10.0 

day 

Ps 

200.0 

lbm/ft° 

r 

1.41 

- 

4^ini  — 

0.12 

- 

We  use  the  relative  density  of  the  gas,  together 
with  the  pressure  and  temperature  values,  to  deter¬ 
mine  the  other  properties  of  the  gas  [8,  20,  26]. 

4.1.  Mesh  refinement 

Table  2  shows  the  numbers  of  cells  used  in  the 
generation  of  the  different  computational  meshes  ap¬ 
plied  in  the  study  of  refinement  of  the  computational 
mesh. 

As  the  production  time  increases,  we  can  see 
(Fig.  4)  that  the  pressure  curves  of  the  well  approach 
each  other.  However,  for  the  initial  instants  of  time, 
there  is  no  overlap  of  values  due  to  the  increased  ef¬ 
fect  of  the  well  known  numeric  artifact  [10,  11]. 
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The  numerical  artifact  appears  as  a  consequence 
of  the  well-reservoir  coupling  technique,  that  assumes 
a  steady-state  flow  regime  near  the  well  [23,  24],  Nev¬ 
ertheless,  its  effects  do  not  compromise  the  results 
for  the  advanced  time  instants.  So,  having  verified 
the  numerical  convergence.  Mesh  5  was  chosen  as 
the  standard,  given  the  significant  reduction  of  the  nu¬ 
merical  artifact  for  times  higher  than  one  day  of  pro¬ 
duction. 


Table  2:  Meshes. 


Mesh 

nx 

Uy 

1 

21 

21 

2 

41 

41 

3 

81 

81 

4 

161 

161 

5 

321 

321 

6 

641 

641 

We  emphasize  that  de  Souza  [9]  validated  the  sim¬ 
ulator  used  in  the  present  work  by  comparing  its  re¬ 
sults  with  those  of  the  commercial  simulator  IMEX  [6] 
for  Model  1 . 


t  (days) 


Fig.  4:  Results  for  mesh  refinement. 

4.2.  Sensitivity  analysis  and  model  comparison 
After  conducting  the  mesh  refinement  study,  we 
performed  a  sensitivity  analysis  considering  the  dif¬ 
ferent  models  for  the  single-phase  flow  of  natural  gas, 
incorporating  (or  not)  the  slippage  and  gas  adsorption 
phenomena.  We  must  remark  that  for  all  the  studied 
cases,  the  Kn  remained  within  the  validity  range  of 
the  slip  flow  regime  [4]. 

Initially,  we  compare  the  results  obtained  with 
Models  1  (Fig.  5)  and  6  (Fig.  6)  in  simulations  by 
varying  the  absolute  permeability.  We  can  see  in 


Fig.  5  that  for  the  highest  values  of  absolute  perme¬ 
ability,  there  are  the  smallest  pressure  variations  in 
the  well.  In  the  graphs,  the  sloping  lines  are  charac¬ 
teristic  of  the  transient  regime  in  a  porous  medium. 
In  Fig.  6,  again,  for  the  highest  permeabilities,  there 
is  the  smallest  pressure  variation  in  the  well.  We 
also  observed  the  occurrence  of  the  transient  regime. 
However,  compared  with  Fig.  5,  it  appears  that  the 
results  for  Model  6  show  a  lower  pressure  drop  than 
those  for  Model  1 ,  for  the  same  absolute  permeabil¬ 
ity,  as  a  result  of  higher  apparent  permeability.  It  also 
leads  to  a  shorter  duration  and  magnitude  of  the  nu¬ 
merical  artifact.  Besides  that,  we  should  note  that  a 
lower  absolute  permeability  leads  to  a  higher  Kn  [8], 
thus  leading  to  a  lower  pressure  drop,  in  contrast  to 
the  effect  of  the  absolute  permeability  in  the  modified 
Darcy’s  law. 


Fig.  5:  Permeability  variation  for  Model  1. 


Fig.  6:  Permeability  variation  for  Model  6. 


About  temperature  variation  for  Model  6,  see 
Fig.  7,  its  increase  leads  to  an  augmentation  in  the 
Knudsen  number  (via  A)  and,  consequently,  to  a 
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growth  in  the  permeability  apparent  (facilitating  the 
flow),  which  contributes  to  a  lower  pressure  drop. 
However,  the  viscosity  of  the  gas  rises  with  tempera¬ 
ture  too.  Thus,  it  also  causes  an  increase  in  pressure 
drop  due  to  the  higher  resistance  to  flow.  Therefore, 
there  is  a  contraposition  of  effects  that  leads  to  dif¬ 
ferent  trends  about  pressure  variation,  with  growth  in 
pressure  drop  for  increasing  temperature  being  preva¬ 
lent  in  the  tested  cases. 

We  can  see  in  Fig.  8  the  influence  that  Vl  has  on 
the  term  source  (due  to  adsorption)  and  the  apparent 
permeability.  We  noticed  that  the  curves  are  close  in 
the  initial  moments  and  that  at  a  later  time,  there  is  a 
gap  between  them.  For  the  higher  volume  of  Lang¬ 
muir,  there  is  a  higher  amount  of  adsorbed  gas  that 
can  be  released,  favoring  the  maintenance  of  produc¬ 
tion  (less  pressure  drop).  We  know  that  the  increase 
of  Vl  causes  a  reduction  in  the  apparent  permeability. 
However,  with  the  release  of  the  adsorbed  gas,  the 
effect  of  Vl  on  k'a  decreases. 


Fig.  7:  Temperature  variation  for  Modei  6. 


Fig.  8:  Variation  of  Vl  for  Modei  6. 


In  the  last  case  studied,  we  compare  the  six  flow 
models,  and  we  show  the  results  in  Fig.  9.  Indeed, 
we  should  point  out  that  we  have  not  found  Model  5, 
ka  =  k'a  and  Tg  ^  0,  in  the  literature.  Nevertheless, 
it  could  be  used,  for  example,  in  reservoir  simulations 
where  the  effects  of  adsorption  are  considerable,  al¬ 
though  we  can  neglect  those  of  slippage. 

When  comparing  Models  1  and  2,  we  observe 
that  the  incorporation  of  the  slippage  effect  leads  to 
a  lower  pressure  drop  due  to  the  higher  values  of  ap¬ 
parent  permeability.  The  curves  start  from  different 
points  and  have  different  slopes,  in  the  region  of  the 
transient  flow,  due  to  the  difference  in  apparent  per¬ 
meabilities.  In  the  specific  case  of  Models  1  and  3, 
the  curves  begin  close  with  each  other,  as  the  appar¬ 
ent  permeabilities  are  the  same  (the  absolute).  Fur¬ 
ther,  we  observe  that  the  slopes  are  similar  in  the  re¬ 
gion  of  the  transient  regime  (except  for  the  effects  of 
non-linearities)  due  to  the  properties  of  the  fluid,  with 
a  lower  pressure  drop  for  the  case  where  adsorption 
supplements  production. 


t  (days) 

Fig.  9:  Comparison  between  different  modeis. 

For  Models  1  and  4,  we  perceive  that  adsorption 
and  slippage  effects  favor  the  maintenance  of  pres¬ 
sure.  Because  of  the  slippage  effect,  we  must  correct 
the  permeability  apparent.  As  a  consequence,  the 
slope  of  these  curves  is  not  similar  to  those  of  Mod¬ 
els  1  and  3.  Indeed,  they  are  closer  to  that  of  Model 
2.  However,  the  overlapping  effects  of  slippage  and 
adsorption  did  not  lead  to  values  of  pressure  higher 
than  those  obtained  in  Model  2.  We  can  explain  this 
behavior  by  the  fact  that  Kn  depends  on  pressure  so 
that  for  higher  pressure  values,  there  is  a  decrease 
in  Kn  and  a  higher  resistance  to  flow.  In  the  case  of 
Model  3,  the  overlapping  effects  lead  to  higher  pres- 
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sure  values  after  about  20  days  of  production. 

As  the  apparent  permeability  in  Model  5  is  lower 
than  the  absolute  permeability,  at  least  at  the  be¬ 
ginning  of  production,  the  numerical  artifact  appears 
more  strongly.  In  the  beginning,  the  production  poten¬ 
tial  is  lower,  but  over  time  the  increase  in  Langmuir 
volume  results  in  a  higher  amount  of  adsorbed  gas 
that  can  be  released,  favoring  the  maintenance  of  pro¬ 
duction  (less  pressure  drop).  Even  so,  this  situation  is 
transient,  and  Model  1  has  a  lower  pressure  drop  over 
long  periods.  A  possible  explanation  is  associated 
with  the  pressure  variation  in  the  reservoir  since  it  has 
a  limited  direct  impact  on  ka  and  because  the  model 
only  considers  its  correction  based  on  Vl-  Therefore, 
the  predominant  effect  becomes  that  of  the  accumu¬ 
lation  term,  and  lower  pressure  drops  than  those  of 
Model  1  would  only  be  possible  for  even  higher  val¬ 
ues  of  Vl- 

Finally,  for  Model  6,  the  difference  in  pressure  drop 
is  highlighted  when  slippage  and  adsorption  phenom¬ 
ena  are  not  neglected,  as  in  the  case  of  Model  1 .  Al¬ 
though  the  effects  are  opposite  concerning  the  vari¬ 
ation  in  apparent  permeability  (predominant  in  the 
cases  studied  here),  adsorption  contributes  to  sus¬ 
taining  production  through  the  accumulation  term. 

V.  CONCLUSION 

As  expected,  we  showed  the  positive  influence  of 
adsorption  once  we  have  more  gas  available  for  pro¬ 
duction  as  a  consequence  of  this  phenomenon,  de¬ 
spite  its  negative  impact  on  the  calculation  of  the  ap¬ 
parent  permeability.  On  the  other  hand,  we  also  re¬ 
marked  the  benefits  of  the  slippage  effect  on  increas¬ 
ing  the  apparent  permeability,  facilitating  flow  through 
the  reservoir.  Besides,  we  were  able  to  capture  the  in¬ 
fluence  of  non-Darcy  models  in  the  well-reservoir  cou¬ 
pling  by  analyzing  the  pressure  variation  of  the  pro¬ 
ducing  well  during  production. 

Therefore,  it  was  clear,  given  the  results  obtained, 
that  for  single-phase  gas  flow  in  shale-type  reservoirs, 
failure  to  consider  both  phenomena  of  slippage  and 
gas  adsorption  may  result  in  incorrect  values  for  pres¬ 
sures  in  the  well  and  the  gas  reservoir.  Such  a  fact 
may  lead  to  mistaken  decision  making.  For  exam¬ 
ple,  we  can  invest  resources  in  the  production  of  a 
deposit  that  will  present  less  pressure  on  the  produc¬ 
ing  well  than  expected  or  failing  to  invest  in  a  reservoir 
that  would  have  favorable  production  conditions.  As 


we  have  seen,  this  can  happen  if  we  do  not  consider 
these  two  effects.  Further,  we  could  perform  simula¬ 
tions  resulting  in  production  flows  closer  to  real  con¬ 
ditions  if  we  take  into  account  both  slippage  and  ad¬ 
sorption  phenomena. 

In  the  case  of  shale  gas  reservoirs,  due  to  their 
characteristics,  we  must  not  forget  that  their  feasibil¬ 
ity  is  possible  through  the  use  of  horizontal  producing 
wells  and  hydraulic  fracturing. 
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